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Cooperative binding affects many processes in biology, but it is commonly addressed only in equi¬ 
librium. In this work we explore dynamical cooperativity in driven systems, where the cooperation 
occurs because some of the bonds change the dynamical response of the system to a regime where 
the other bonds become active. To investigate such cooperativity we study the frictional binding 
between two flow driven surfaces that interact through a large population of activated bonds. In 
particular, we study systems where each bond can have two different modes: one mode corresponds 
to a fast forming yet weak bond, and the other is a strong yet slow forming bond. We find consid¬ 
erable cooperativity between both types of bonds. Under some conditions the system behaves as if 
there were only one binding mode, corresponding to a strong and fast forming bond. Our results 
may have important implications on the friction and adhesion between sliding surfaces containing 
complementary binding motifs, such as in the case of cells binding to the vessel walls under strong 
flowing conditions. 


Many biological functions depend on the ability to ar¬ 
rest objects traveling in the blood stream by binding 
them to the vessel wall. Examples include the binding 
of platelets to stem hemorrhage [1, 2] and the binding 
of immune cells to the vessel wall prior to their entry 
into adjacent tissue [3, 4]. These bonds are formed be¬ 
tween dedicated proteins existing in the blood stream, on 
the platelets or immune cells, or on the vessel wall. In 
order to arrest the flowing cells, several ligand-receptor 
bonds need to form, particularly at high flow rates. Such 
ligand-receptor bonds, as for example selectin and its lig¬ 
and in leukocytes [3], or glycoprotein Gpiba in platelets 
and the A1 domain in Von Willebrand Factor (VWF) 
[2, 5, 6], are protein complexes that contain elementary 
bonding units such as salt bridges and hydrogen bonds, 
as well as hydrophobic interactions. To simultaneously 
activate multiple elementary bonds or effective hydropho¬ 
bic interactions requires configurational registry between 
the protein and ligand, which means that the strongest 
bonds tend also to be the slowest. 

Interestingly, recent developments have shown exceed¬ 
ingly rapid adhesion of platelets to VWF; the time 
scales are the fastest that have been reported for ligand- 
receptor interactions [7-9]. This behavior is still far from 
understood, yet it is critical for the formation of plugs 
[2, 10-12]. Here, we demonstrate that a system in which 
we combine multiple types of bonds can exhibit such be¬ 
havior. In particular, we study a simple model of how an 
object embedded in a fluid and flowing over a substrate 
arrests under the action of a cascade of binding events. 
We are interested in a scenario in which a single kind 
of bond is not sufficient to arrest the system. The first 
kind of bond is not strong enough, the second is not fast 
enough. In this scenario we study the effect of dynam¬ 
ical cooperativity between the bonds when both types 
are present, and find the conditions under which they 


together confer their favorable properties (being fast or 
being strong) on the system. Our results indicate that 
depending on the degeneracy of the bonds and the dy¬ 
namical properties of each, it is possible to create a sys¬ 
tem which appears as if it contains a single “hybrid bond” 
type that is fast forming yet strong. This is quite rele¬ 
vant in fast sliding conditions such as blood flow, where 
common biological bonds may be too slow to significantly 
affect the dynamics. Such a scenario is particularly im¬ 
portant during uncontrolled formation of platelet-VWF 
aggregates in stenotic regions that can lead to stroke or 
heart infarcts [9, 13, 14] . 

Our work sheds light on the importance of weak bonds 
in out-of-equilibrium conditions. The existence of weak 
bonds has long been speculated to have very important 
implications in binding because they may increase the 
apparent on-rate substantially. For example, it was re¬ 
cently pointed out that a variety of binding sites can be 
present on protein surfaces, and that even weaker bind¬ 
ing regions can significantly affect the kinetics of binding 
[15-18]. 

Our model is related to the method called adhesive 
dynamics [19-23], but we have simplified the flow prob¬ 
lem in order to focus solely on the interplay between the 
bond types, which to the best of our knowledge has not 
been addressed before. The model is sketched in Fig. 1. 
It consists of a substrate and a rigid object immersed in 
the fluid flow above the substrate, and the bonds between 
them. The substrate can for example represent the wall 
of a blood vessel and the rigid object a platelet in the 
blood stream. The bonds between the substrate and the 
rigid object are modeled as springs. The intact bonds 
contribute a friction force that resists the motion of the 
object, slowing it down or arresting it depending on the 
bond properties and flow conditions. The equation of 
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FIG. 1. (Color online) Sketch of the model system. A rigid 
object interacts with the substrate though the formation and 
breaking of bonds. Each bond can be in a broken, a strong 
or a number of weak configurations. 

motion of the rigid object is 

MX = riiVfiuid -X)+ Fb, (1) 

where M is the object’s mass, X is its acceleration (the 
dots denote time derivatives), rj is the effective viscosity 
of the object in the fluid, 14 uid is the speed of the fluid 
with respect to the substrate, X is the object’s instanta¬ 
neous speed with respect to the substrate and Fi, is the 
force from the intact bonds, 

N 

Fb = J2ft^- ( 2 ) 

N is the number of bonds that can form. Broken bonds 
contribute zero force. The force in an intact bond is 
fi = K{li — Iq), where k is the bond stiffness, h is the 
length of the bond and Iq is an equilibrium bond length; 
h = y/xf + where Xi is the elongation of the bond 
along the x direction and h is the vertical separation of 
the object and the substrate. The projection of fi along 
the direction of the fluid flow is = —fiXi/k. When a 
bond is intact the elongation of the bond occurs with the 
same speed as the motion of the object, ii = X. When 
the bond breaks, its elongation returns to zero, controlled 
by the parameter A as follows 

._iX , intact, 

* \ —Xxi , broken. 

The rupture and formation of bonds is controlled by 
two rates that we label and A:°". These depend on the 
type of bond being broken or formed (weak or strong). 
In each time step At of our simulation an intact bond 
breaks with probability k°^At. We draw a random num¬ 
ber from the uniform distribution f G [0, 1] and break 
the bond if < k°^At. For bond formation, we first 
draw a random number f to determine which type of 
bond to attempt to form, with each bond type weighted 
by their relative abundance. Thus, if there are 5 differ¬ 
ent weak bonds and 1 strong bond, Atype = [5, 1], the 


attempt is for a weak bond if ^ < 5/(5 -I- 1) and for a 
strong bond otherwise. This can be generalized to any 
number of bond types. Then, another ^ is drawn and 
the bond forms if ^ < Atypesfc™At. Multiplying by the 
sum of bond types, Atypes = X]i-^type,i) ensures that 
the bond formation attempt probabilities are additive, 
so that the probability to form each bond is independent 
of the number of bond types being considered. The rates 
k°’^ are thus characteristic of each type of bond, while 
the number of potential bonds of each type is encoded in 
Atype. Even though bond formation attempts are addi¬ 
tive, however, each bond can only be of one type at any 
given time (or the bond can be broken). 

For the bond rupture rates we use the simple Bell-type 
expression [24] 

= /eg® exp (/3/i Ax), (4) 

where fcp® is the bond rupture rate in the absence of any 
force on the bond, /3 = 1/(/cbT) and Ax is the length 
to the transition state from the bottom of the well. For 
biological protein-protein complexes this scale is on the 
order of a nanometer. Filippov et al. [25], studying a 
similar model, introduce another rupture rate expression 
for what they call strong bonds, but we have used equa¬ 
tion (4) for both types of bond, assigning different /eg® to 
weak and strong bonds. 

The rate of bond formation when the object is station¬ 
ary with respect to the substrate is for each type of 
bond. The rate of bond formation when the object is in 
motion is the product of the stationary rate and a func¬ 
tion g{) that reduces the rate. This function effectively 
introduces a limiting timescale for the formation of such 
bonds. This is natural in all systems, as there is always a 
relaxation timescale upon which certain binding configu¬ 
rations can be reached. One can consider this timescale 
as an effective attempt rate. In the case of proteins, 
one expects such a time scale to be in the millisecond 
range. Interestingly, in the case of the interaction be¬ 
tween Gpiba and Al such interaction appears to occur 
within microseconds. When two surfaces slide past each 
other, the time available to form the bond is r = a/|A|, 
where a is the lateral extent of the reaction site. We 
assign a characteristic time tq (or equivalently a charac¬ 
teristic velocity Vg) for each bond type, tq = a/Vo. The 
full rate of bond formation is 

, (5) 

with At a parameter that controls the sharpness of the 
reduction in with |A|. We have used g( “ 

tanh(^^^^;^). We express length, time, and force in units 
of h, M/rj, and ksT/h, respectively. 

We set as initial conditions that the object is moving 
with the fluid with no intact bonds, X = Efluid, h = h'^i- 
We then integrate the motion of the rigid object forwards 
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FIG. 2. (Color online) Example time evolution of object speed 
and bond properties when (a) only weak (but fast) bonds are 
possible, and when (b) both weak (but fast) and strong (but 
slow) bonds are possible. Top panels: object speed and aver¬ 
age force in bonds. Bottom panels: fractions v of unformed 
(-), weak ( ) and strong (-) bonds. In this paper N = 300, 
r; = N/A:, M = N, k = [10, 10] (for [weak, strong] bonds), 
lo = h = I, \ = O.dry/M, kf = [0.2, 0.001], A:g“ = [0.1, 0.1], 
Eo = [K^“^ l/o“] = ]1, 0.1], a = [1, 1], Ar = [0.1, 0.1], 
/3 = 1, Aa; = [1, 1], At = 0.01. In this figure Ntype = [5, 0] 
(a) or iVtype = [5, 1] (b), 14uid = 0.8. 


in time using the leap-frog method with time step At. 
The resulting motion has three distinct regimes: (i) the 
object moves with the flow with very little friction from 
the substrate and the average speed of the object {X) « 
14uid, (ii) significant friction from the substrate slows the 
object down to {X) markedly smaller than Vnuid, and 
(hi) the object comes nearly to rest. The regime selected 
depends on the properties of the bonds available and on 
the flow speed Tsuid- 

Figure 2a shows a typical evolution for weak bonds. 
The population of bonds quickly stabilizes around an 
equilibrium value, with bonds being broken and formed 
continuously. The object slows down to an average speed 
{X) « 0.1, markedly smaller than Fjiuid = 0.8, but far 
from coming to rest. Figure 2b shows a typical evolution 
when a strong bond option is added; the settings are oth¬ 
erwise the same as in Fig. 2a. The beginning of the dy¬ 
namics is the same. The weak bonds form and the speed 
of the object is reduced. Notice that for the first 18 or so 
units of time, the number of strong bonds remains neg¬ 
ligible. Then, as the object’s speed falls to A « 
and the probability to form strong bonds increases sig¬ 
nificantly, the strong bonds take over and the object’s 
speed eventually falls to A « 0. For additional data on 
steady-state bond fractions see Supplemental Material at 
the end of this file. 

Figure 3 shows the average object speed (A) as a func¬ 
tion of fluid flow speed 14uid for a range of scenarios for 
what bonds are available. Figure 3a shows three cases. 


(i) If only strong (slow) bonds exist, the object moves 
as fast as the fluid for all but the lowest fluid speeds. 
For Vfluid < strong bonds are formed and bring 

the object to rest, (ii) If only weak (fast) bonds ex¬ 
ist, the object is slowed down by the weak bonds for 
f4uid ^ but the bonds are unable to bring the ob¬ 

ject to rest, (hi) If weak (fast) and strong (slow) bonds 
both exist, their combined effect brings the object to rest 
for Ffluid ^ Figure 3b shows the combination of 

I, 3, 5, or 7 weak bonds with 0 or 1 strong bond. When 
there are fewer weak bonds they are less effective at slow¬ 
ing the object down, and the transition to arrest in the 
weak plus strong case moves to lower values of Ffluid- 
Note that for Atype = [7, I] the combined dynamics is 
equivalent to what we would expect from having bonds 
that are strong and fast. Thus, we can say that the fast 
bonds lend their speed to the slow bonds, or equivalently 
that the strong bonds lend their strength to the weak 
bonds. This synergistic effect is what we termed dynam¬ 
ical cooperativity. 

For the conditions that lead to the object coming to 
rest. Fig. 4e shows the arrest time as a function of t4uid- 
In all our simulation the arrest occurs after the fraction 
of strong bonds starts growing. We use the first time 
when r'strong bonds > 0.2 as an estimate for the arrest 
time, because this is numerically more robust than a test 
on A, which can fluctuate around or near zero, see Fig. 2. 
For each combination of available bonds there is a tran¬ 
sition speed i/^™sition above which strong bonds do not 
form. The arrest time diverges as I4uid —t 
from below. Away from the divergence, there is only 
a weak dependence of the arrest time on kjiuid- This 
can be understood from Fig. 3b. When the weak bonds 
alone lead to an average speed (A) < the object 

first slows down under the influence of the weak bonds, 
then the strong bonds take over and the object arrests. 
This process depends only weakly on I4uid as long as 
bfluid < V]]"®®'’'. When the weak bonds alone lead to 
(A) ^ however, the strong bonds only form when 

fluctuations in the object’s speed momentarily bring A 
below Thus the arrest time becomes dependent 

on the extreme value statistics of the fluctuations in A. 
Figure 4a-d show the distribution of arrest times for se¬ 
lected conditions. Away from j^ransition (a^^j) arrest 
times are well fitted by normal distributions. Closer to 
(bjC) the distributions of arrest times change 
and become strongly asymmetric. 

In this work we have demonstrated that the kinetics 
of binding is strongly modified when one has dynamical 
cooperativity in which two different types of bonds ex¬ 
change their dynamical properties. In our case we showed 
that under some conditions one can obtain complete “hy¬ 
bridization” of the properties, and the system appears 
as having bonds that are fast and strong, although one 
starts from a mixture of slow and strong bonds together 
with weak and fast bonds. Such cooperativity is interest- 
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FIG. 3. (Color online) Average steady-state speed of rigid ob¬ 
ject vs fluid flow speed Vkuid for systems with different bonds 
available, (a) Strong only, weak only, and a combination, 
(b) Weak only and the corresponding weak plus one strong 
bonds. The gray dashed line shows X = Vfluid as a guide to 
the eye. Legends show values of Atype- {X) was measured 
over the second half of simulations lasting for 150 units of 
time. Points near the transition from arrest to non-arrest are 
averages over ca. 20 simulations. 


ing because it relies heavily on the degeneracy of the weak 
bonds, and not on the strong bonds. In other words, it 
is dominated by the kinetics of the bonds, and not the 
equilibrium properties. 

In conclusion, for strongly out-of-equilibrium situa¬ 
tions we have found that the presence of the weak bonds 
is crucial, and enables the system to achieve states (such 
as the arrested state) that would otherwise be unattain¬ 
able. This leads us to suggest that experimental obser¬ 
vations of apparently fast and strong bonds may in fact 
be due to dynamical cooperativity rather than to a single 
bond type. The question is how to detect the presence 
of the weak bonds. In principle, this could be possi¬ 
ble by studying the rolling of beads under an applied 
torque as a function of the rotational frequency [26]. If 
the timescales are well separated, one should see two dif¬ 
ferent regimes of friction dominated by the different bond 
types. Clearly, such dynamical cooperativity would be of 
much importance in situations where one would like to 
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FIG. 4. (Golor online) Estimated object arrest time, (a-d) 
Distributions of arrest times for selected simulation settings 
(arrows originate on the Atype = [3, 1] and [5, 1] lines). The 
fits in (a) and (d) are normal distributions. Each histogram 
is based on ca. 2000 simulations, (e) Estimated arrest time 
vs Ffluid- Golors and markers are the same as in Fig. 3. For 
the points close to where each line diverges, the total sim¬ 
ulation time was increased to 300 units of time. Error bars 
show one standard deviation away from the mean. Away from 
the error bars are smaller than the markers. Each 
point is based on 40 simulations. 

“recycle” different binding motifs, and not synthesize de 
novo a special binding partner with the prescribed ki¬ 
netic attributes. Thus, we expect this behavior to be 
important in many diverse areas of lubrication, friction 
and adhesion. Furthermore, confirming and characteriz¬ 
ing such weak interactions could lead to progress in our 
understanding of the equilibrium binding properties of 
proteins [15-18]. 
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SUPPLEMENTAL FIGURE 


Figure S5 shows the average bond fractions in the 
steady state reached at the end of our simulations. Each 
panel combines Vype = [^weak, 0] and [iVweak, 1], with 
-^weak = 1, 3, 5. We measured the average bond fractions 
in the last 10 units of simulation time for the same simu¬ 
lations as in Fig. 3. Where there are multiple simulations 
with identical settings, the average over simulations was 
taken as well. 


FIG. S5. Average bond fractions reached in steady state. In 
the simulations with only weak bonds, the average fraction of 
formed bonds has an appreciable negative slope thronghout 
the range where bonds are formed. In the simulations with 
strong bonds as well, the fraction of bonds formed stays nearly 
constant up to s-nd then drops sharply. This is 

consistent with our observation in the main text that strong 
bonds arrest the object, but weak bonds do not, because the 
bond fraction depends on the slider speed, which becomes 
independent of Vkuid as long arrest occurs. Legend: fractions 
of unformed (-), weak ( ) and strong (-) bonds for Wype = 
[1, 0] (b), [1, 1] W, [3, 0] (o), [3, 1] (.), [5, 0] (□) and [5, 1] (.). 


SUPPLEMENTAL CODE 


Our computer code is available from arXiv as a sepa¬ 
rate download. 














